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■ Abstract 

We present the complete set of analytical solutions of the geodesic equation in the five- dimensional Myers- 
. Perry space-time with equal rotation parameter in terms of the WeierstraB' elliptic and WeierstraB' zeta and 

' sigma functions. We study the underlying polynomials in the polar and radial equations which depend on the 

parameters of the metric and conserved quantities of a test particle and characterize the motion by their zeros. 

We exemplify the efficiency of the analytical method on the orbits of test particles. 

fcnJ 1 Introduction 

ox) 

The Myers-Perry [1, space-time is a higher dimensional generalization of the Kerr space-time which describes rotating 
black holes in more than four dimensions. Interest in higher dimensional gravity is a result of the search for the 
theory of quantum gravity which should describe the Universe at its beginning. With (Super) String and M- theory 
as candidates, it finds applications in the AdS/CFT correspondence. Another interesting question is the creation of 
black holes at the LHC which from an energetical point of view is only in higher dimensions possible (associated with 
£C) ' a change of the fundamental value of the Planck mass in higher dimensions). A review of the higher dimensional 
black hole solutions in vacuum and in supergravity theory is presented in [2[ |3] and a discussion of the physical 
properties of higher dimensional black holes is addressed in [4 . General Kerr-de Sitter and Kerr-NUT- AdS metrics 
in all dimensions are presented in [SJ [6] . A very good motivation to study the higher dimensional black holes is 
proposed in [7]. 

In the present work we study the motion of test particles in the 5D Myers-Perry space-times. For simplicity of 
the analysis we have chosen equal values of the rotation parameters. The case with unequal rotation parameters 
is currently under study [8] . The motion in the Kerr space-time is presented in [9j [10] . Some special cases of the 
motion in 5D Myers-Perry space-times are studied in [IT], where it was shown that there are no stable circular 
orbits in the equatorial place, and in [12] scattering and capture of particles by 5 dimensional rotating black holes 
are investigated. Geodesic stability of circular orbits in singly-spinning d-dimensional Myers-Perry space-times is 
studied in [13] . Properties of the motion of charged test particles around a black hole immersed into a magnetic field 
are studied in [14 . Separability of the Hamilton- Jacobi equations of motion in the Kerr-de Sitter and Myers-Perry 
space-times in all dimensions are investigated in [15j [16] and in [17j [18] [19] [20] [21] . Analytical solutions of the 
geodesic equations in the Kerr-de Sitter space-time in 4D are studied in [22] [23] . A general analytical solution of 
the geodesic equations in terms of the hyperelliptic theta- and sigma- functions in the Myers-Perry space-time with 
one rotation parameter is presented in [24] . 

In this paper after the introduction of the equations of motion in Sec. [2] we discuss in Sec. [3] the properties of 
test particle motion on the basis of the polar and radial equations and derive in Sec. 2] the complete set of analytic 
solutions of the geodesic equation in the Myers-Perry space-time in 5 dimensions for equal values of the rotation 
parameters in terms of the WeierstraB' elliptic, zeta- and sigma- functions. In Sections [5] and [6] we plot the orbits 
for chosen sets of parameter values and discuss observable quantities. 
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2 The geodesic equation 

'ic of a five-dimensional rotating 1 

f dlJ 2 \ V 2 9 

ds 2 = -dt 2 + p 2 f — + d$ 2 J + a sin 2 tidip 2 + /? cos 2 # 2 + -± (dt + a sin 2 tidtp + b cos 2 #) , (1) 



The metric of a five-dimensional rotating black hole in Boyer-Lindquist coordinates has the form [1, 5 , 25 , 11 

,2 \ ^2 



where p 2 = u + o 2 cos 2 $ + 6 2 sin 2 $ and a = + a 2 , f3 = u ^ b 2 ^ A = a/3 — TqW. Here we have introduced the 
coordinate u such that: u = r 2 and r is a radial coordinate. The parameter r 2 , is proportional to the mass of the 
black hole. 

The metric is singular at A = and p 2 = 0. There are also angular coodinate singularities at $ = and $ = ^. 
Calculating the Kretschmann scalar one sees that the surface p 2 = is a scalar curvature singularity [25] since the 
Kretschmann scalar diverges at p 2 = 0. Thus, for a = b the singularity is located at u — —a 2 [25] as shown in FigJTJ 
The zeros of A define the horizons which for a = b yield 



u± 



rl - 2a 2 ± r 



4a 2 



For 4a 2 = Tq the horizons merge and the space- time becomes extreme. 
The ergosphere defined by g u = is located at 



^stat =rl - a 2 . 



(2) 



(3) 



and does not depend on the angle $ for a = b. 
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Figure 1: Position of the singularity and the horizons in the Myers-Perry black hole with equal rotation parameter a — b. 



For test particles or light (g^ v u^u v = —5, S = 1 for massive test particles and S — for light) moving on a 
geodesic the energy E and angular momenta and ^ are conserved. 

It is straightforward to see that in the Myers-Perry space-time in 5 dimensions the Hamilton-Jacobi equation 



2 ds = as_as_ 

dr dx^ dx v 



(4) 



where r is an affine parameter along the geodesic, is separable and yields for each coordinate a corresponding 
differential equation 
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Here we used £ = Eaj3 + a<&j3 + b^a and introduced the Mino time 7 through p 2 dj = dr [26]. We also defined 

£ 2 
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A((E — 5)u — K + V) 



] af3 



6 = (E 2 - 6)(a 2 cos 2 + b 2 sin 2 tf) + K 



$2 



* 2 



sin 2 $ cos 2 1? ' 



(10) 
(11) 



where V = (a 2 — b 2 ) — 7r)' ^^ ie se P ara ti° n constant K is known as Carter constant. 
In the following chapters we consider the case of equal rotation parameters: a = b. 



3 Complete classification of geodesies 

Consider the Hamilton- Jacobi equations ((5])- ((9]). The properties of the orbits are given by the polynomial U (fTUj) 
and the function ([TT]) . The constants of motion (energy, angular momenta and separation constant) as well as 
the parameters of the metric (mass, rotation parameter) characterize these polynomials and, as a consequence, the 
types of orbits. In this section we discuss the motion in Myers-Perry in terms of the properties of the underlying 
polynomial U and the function 0. 
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Table 1: Types of orbits in the Myers-Perry space-time for the rotation parameter a = b. The horizontal bold lines represent the orbits and the 
bold points - the turning points. The thin vertical double dashed line represents the horizons and the thick vertical dashed line the 
singularity at u = —a 2 . In special cases the turning points lie at the horizons (±), or lead into the singularity (s). The regions (l)-(4) 
are related to the E-<& diagrams for the u-motion in Fig. [3] 



3.1 The ^-motion 

We rewrite the function O (jlip in the following form: 

t-Hi ^ i „ tx. T x2 ($cos 2 tf + *sin 2 tf) 2 

cos^ v sin v 

9 _ _ , x tx9 ($cos 2 ^- ^sin 2 ^) 2 , . 

= (E 2 -5)a 2 +K-{$ + ^) 2 - K - 2 q • 2 n ~ • ( 12 ) 

COS z 17 Sin 17 

In order to obtain from Eq. (j6j) real values of the coordinate d we have to require 6 > 0. This implies 

ci :=(E 2 -5)a 2 +K-(<S>-^) 2 >0 

c 2 := (E 2 - S)a 2 + K - ($ + ^) 2 > . 1 ^ 

With the new variable £ := cos 2 Eq. ((6]) turns into the equation 

^f) =6^ with 6^:=a^ 2 + ^ + c^, (14) 



with a simple polynomial of second order on the right hand side, where 



^ = -(c 1 + ($-*) 2 ) = -(c 2 + ($ + *) 2 ) , (15) 

^ = K + a 2 (E 2 -5)+^ 2 -$ 2 , (16) 

| = "* 2 - (17) 

Since c\ > we have < 0. This means that can be positive if and only if there are real zeros of 0£. The 
polynomial plays the role of an effective potential for the ^-motion. The real and positive zeros of define 
the angles of two cones which confine the motion of the test particles. (A similar feature appears in Taub-NUT, 
Kerr-(de Sitter) and Reissner- Nordstrom (for charged particles) space-times [27 ] 123 ] [28].) 

The discriminant D = b 2 — 4a^c^ of the polynomial Be can be written as D = 16ciC2- The existence of real 
zeros of requires D > 0. And with the conditions ([T3]) this is fulfilled. The inequalities (fT3|) impose limitations 
on the parameters E, <l>, ^ and K for some given a. 

The polynomial 0^(£) in ([H]) has roots £i 5 2 = — ^ (b% ± y/D^j, and describes a parabola with the maximum 
a ^ ( — — ia^)* """ n s P ec ^ cases the roots £i 5 2 are: 

1. ci = ^1,2 = ^z$. For $ = v[/ the orbit lies at = f . 

2. c 2 = ^1,2 = r^L-. For $ = ^ the orbit lies at = f . 

3. ci = C2 = => 2 cases are possible: £1,2 = 1 for $ = and £i ? 2 = for \£ = 0. Thus, the motion is in the 
r d = plane or r d = ^ plane. 

Fig. [2] shows the dependence on the parameters a, \I/ from the function ([T2]) . Only in the region (d) 
ci, C2 and D are positive, in all other regions either D is negative or one of the conditions ([T3]) or both is/are not 
fulfilled. On the lines the discriminant D vanishes and the motion is characterized by a constant angle For small 
values of the rotation parameter a the curves are less bent and for a — the curves become straight lines. For 
growing values of the separation constant K the distance between the regions (b) grows. By varying ^ we change 
the size of the region (d). 

3.2 The ^—motion 

The right hand side of the differential equation (j5j) has the form U = Y^=o ^°^ u% w ^ n ^ ne coefficients bi for a = b: 

b 3 = 4 (E 2 - 5) 

b 2 = Sa 2 (E 2 - 5) + 4 ( 8rl - K) 

(18) 

h = 4a 4 (E 2 - 5) + AK{rl - 2a 2 ) + Srl(a 2 E 2 + aE($ + *)) 
6 = 4a 2 (rl (aE + ^ + ^) 2 - Ka 2 ) . 

The list of the possible orbits reads (For an explanation of the names of the orbits see [29]): 

1. Escape Orbits (EO) in the region [m, oo) with 1*1 > i*+. The EO does not cross the horizons. 

2. Two-world escape orbits (TEO) in the region [1*1,00) with \m\ < u- 

(a) Terminating escape orbits (TEO s ) in the region [— a 2 ,oo). The test particle comes from infinity and 
disappears into the singularity (u = —a 2 ). 

3. Periodic bound orbits (BO) in the region u £ [1*1,1*2] with 1*1 > 1*2 and 1*1,1*2 < 1*- (planetary orbits). 

(a) Terminating bound orbit (TBO) in the region u £ [—a 2 , 1*1] with 1*1 < i*_. The turning point of a bound 
orbit coincides with the singularity for some energy value. 




(d) a = 0.4, K = 3, ^ = 0. (e) a = 0.4, K = 3, ^ = 1.5. (f) a = 0.4, K = 10, ^ = 0.5. 



Figure 2: The diagram for the ^-motion for different values of a, K, and ^. In the region (b) the discriminant is negative D < 0. In the 
regions (a), (d) and in some cases in the region (c) the discriminant is positive > 0, but only in the region (d) the conditions (|13J) 
are fulfilled. On the lines the discriminant vanishes. Moreover, on the dashed lines one of the quantities ci,2 is negative and the 
motion is not allowed. 



4. Many-world periodic bound orbits (MBO) in the region u £ [^1,^2] with u\ < U- and U2 > u+. 

(a) Also for the MBO it is possible that the turning point concides with the singularity and the orbit turns 
to a TBO in the region u £ [— a 2 ,^i], with u\ > u + . 



Analysing the polynomial U in the equation ([5]) we construct the E-Q diagrams in Fig. [3] for the ^-motion. The 
diagrams illustrate how the change of the values of the parameters a, K, and \I/ influences the number of real 
zeros of the polynomial U . From Table [T] we read the possible types of orbits in these regions. In the light grey 
region are at most 3 real zeros possible. If all three roots belong to the interval [—a 2 , 00), then according to the 
Table[T] we have regions (2) or region (4) depending on the values of the parameters. It can happen that only 2 of 
3 zeros lie in the interval [—a 2 , 00). This corresponds to the region (1). There is only 1 real zero in the region (3) 
painted in dark grey. The dashed area is forbidden by the inequalities ([13j) . 

Region (1): Two real zeros in [—a 2 , 00). According to Table [T] possible orbits are MBO. Special cases are B± when 
the turning points lie at the horizons, B_ and B + when one turning point lies at the inner or outer 
horizon correspondingly. A very specific orbit has the left turning point merged with the singularity 
(B s ). 

Region (2): Three real zeros in [—a 2 , 00). Possible orbits are MBO und EO with special cases C+, C_, C s and C s+ . 

Region (3): One real zero in [—a 2 , 00). Possible orbits are TEO for test particles coming from u = +00. The turning 
point may lie at the inner horizon (A_), or merge with the singularity (A s ). 

Region (4): Three real zeros in [—a 2 , 00). Possible orbits are BO and TEO. BO lies behind the inner horizon. The 
turning point of the TEO may lie at the inner horizon (D_), or the left turning point of the BO may 
merge with the singularity (D s ). 



(a) r = 1, a = 0.2, K = 3, ^ = 0.5. (b) r = 1, a = 0.4, K = 3, ^ = 0.5. 
Definition of the regions (1), (2) 
and (3) (see Table [fj 



(c) r = l,a = 0.4, K = 5, * = 
-0.5. 




Figure 3: The diagrams for the u-motion for different values of the parameters a, K, <3? and \&. The light grey region corresponds to three 

real zeros of the polynomial U and the dark grey denotes the regions with 2 real zeros. When plotting an orbit one has to check in 
addition which of these real zeros lies in [—a 2 , oo). The dashed region is forbidden by the inequalities (|13J> . 



Another useful method to characterise the motion of a test particle is to plot the effective potential V G s. We 
define it in the following way 



/ du x 2 



\dh( 

where 



4 (An + a 2 r 2 ) (E - V+) (E - V~ s ) , (19) 



Vt = 

eff Au 



^—^S (atari (* + *) ± \/ A ( K + 5u ) ( Au + a2r o) ~ &a 2 r 2 u (i + *) 2 ) . (20) 



The potentials are shown in Fig.2]for massive and in Fig.[5]for massless test particles. For light the same types of 
orbits as for the massive test particles exist. The motion is forbidden in the grey region. At infinity lim = 

is true. The potentials and meet at the horizons where A = 0: 



a 



^) = --^(* + *), (21) 

a 



V&«-) = -^(* + ¥) 

If $ + ^ = then the ordinates of both points where the potentials meet lie at the i^-axis and the potentials are 
symmetric wrt the n-axis: = -V^ (Fig. |4(a) ). The potentials are also symmetric wrt the change of sign of the 



parameters of the test particle: V^g (£", 3>, \£) = (-£", — — 

There is in general a potential barrier which prohibits a test particle from falling into the singularity (Figs. | ^a)f 
(d)[) . The potential plots are truncated at the singularity. If a part of a potential is truncated at the singularity, the 



remaining part containing no potential barrier is not allowed by the inequalities (fl"3|) as in the plots (e) and|(f)[ But 



as we know from the discussion above it is possible that a turning point of an orbit coincides with the singularity 
which implies that the orbit ends in the singularity. 

Setting u = —a 2 in the polynomial U with the coefficients ([18]) one sees that only for C2 = the singularity can 
be reached: U(u = —a 2 ) = — 4a 2 rgC2. Thus, if we fix \I/ and <£> then only particles with energy 

E C2 = ±^a 2 5-K + ($ + ^) 2 (22) 

can fall into singularity. The orbits of these particles are denoted by the index s in the Table [TJ 

It is easy to see that only in this case falling into the singularity is possible. The potential is truncated at the 
singularity u = — a 2 . At this point it takes the value 

V^{-a 2 ) = ±-Ja^-K + ^ + ^f (23) 

which coincides with ([22]) . Since C2 > from the inequalities ([13]) , the region between E C2+ = V^(— a 2 ) and 

E C2 _ = V e ^(— a 2 ) is not allowed. But the energies E C2 corresponding to the boundary of this region where C2 = 
are allowed. 




Figure 4: Effective potential for massive test particles (5 = 1). The potentials V^ f (red) and (green) couple at the horizons (dashed 

lines). In the grey regions no motion is possible. The dashed region is forbidden from the conditions on the ^-motion defined by the 
inequalities (1131) . The bold black vertical line at u — —a 2 defines the singularity. The dashed lines indicate the horizons U+ and n_. 



4 Solution of the geodesic equation 

Now we present the analytical solutions of the differential equations ©-(O. 




(d) r = l,a = 0.4, K = 10, $ = (e) r = 1, a = 0.4, K = 10, <£ = 
-4, ^ = -0.5. —4,^ = - 0.5. Detailed repre- 

sentation of |(d)| BO. 



Figure 5: Effective potential for massless test particles (5 = 0). The dotted lines in the potential are plotted for comparison and correspond to 
the test particle with non zero mass. One can see that the types of orbits from Table 1X1 are conserved for light. 

The potentials (red) and (green) couple at the horizons (dashed lines). In the grey regions no motion is possible. The dashed 
region is forbidden from the conditions on the ^-motion defined by the inequalities (|13j) . The bold black vertical line at u = —a 2 
defines the singularity. The dashed lines indicate the horizons u+ and u_. 

4.1 Solution of the equation 

The solution of the differential equation (j6j) is given by the elementary function 



0(7) =arccos(±y^- (Vcsin (- V /= ^7 + 7&) - k)) > (24) 

where < and D > are defined in ([14|h 7^ = ^—a^y in + arcsin ^ ^|^ + — ^ , 7i n is the initial value of 7 and 
^i n is the initial value of £. 

4.2 Solution of the ^—equation 

A substitution u = ^ (4y — y) transforms the polynomial U into the standard Weierstraf] form 
where 



V ^J = ^ - 92V - 93 := Ps(y) , (25) 

52 -12"~' 53 = ^8~~76~ _ 216 " (26) 
The differential equation (|25j) is of elliptic type and is solved by the Weierstraf] ' p-function [30] 

y(i) = p (7 -7^; 92, 93) , (27) 



where j f [n = 7i n + ^ = dy — with y- ln = ^fu- in + f§. Then the solution of (j5]) acquires the form 



u = (^P(l-7L\92,93)- j) • (28) 



4.3 Solution of the 9?— equation 

Eq. ([7]) can be simplified by using ([5]) and ([6]) and by performing the substitution £ = cos 2 # and u = ^- [Ay 



3 , 



where 



where 9^ is given in (fl4|) for a = 6, and 



= dy?£ + d(p u , (29) 



M = 7^ 7 ( 3 °) 



d^ = (si) 



where the partial fractions decomposition procedure was applied. Here Ki, i = 1,2 are constants which arise from 
the partial fractions decomposition. These depend on the parameters of the metric and the test particle, and p\ 
and P2 are the zeros of A = A (u = ^ (4y — y)) • 
The integration of ([29]) yields: 

<P - ^in = <P£ + , (32) 

where <^i n is the initial value of the angle p. 



Consider first the dp u . After the substitution y = p(v) with p'(v) = ^/4p 3 (v) — g2p(v) — #3 and integration of 
the differential ([3T]) one gets: 

f v 2 K 

<p u = V — ^ rdv . (33) 



Here v = ^(7) = 7 - 7^ and v in = v(7 in ). And = pj , j = 1, 2. 

The solution of ([33]) is given in terms of the elliptic a and £ functions [28j [27] : 

^ = E-?rS( ^)(« " ^ + ln aJ ri^ ~ ln t 5 ^) > ( 34 > 
^ p '(« PJ ) V ^ + ) CT ^ in + w ^ ) / 

The integral ip^ can be easily found and the solution for < and D > is given by 

, £(7)= ' arctan ^iS&±M^, (ss) 

£in 



2 4$V=a£^/l - «| 



2a 5 £ + 6 ? 



= • (^) 



where 



4.4 Solution of the ^-equation 

Like for the coordinate tp the Eq. dSj) can be simplified by using ([5]) and ([6]) and by performing the substitution 
£ = cos 2 and tx = ^ (4?/ - ^) 

= dijj^ + <i^ n , (37) 

where 



where is given in ([T4]) for a = b. And 



dV» c = — |=j (38) 



^ = ^=E^> (39) 



where analogously the partial fractions decomposition procedure was applied. Here N{, i — 1,2, are constants which 
arise from the partial fractions decomposition. These depend on the parameters of the metric and the test particle. 
And pi and P2 are the zeros of A = A (u = [4y — -y)^. 
The integration of ([37]) yields: 

%j) - ^ in = ^ + lj) u , (40) 

where ^i n is the initial value of the angle ip. 



Consider first the dtp u . After the substitution y = p(v) with p'(v) = ^/Ap 3 (v) — g2p(v) — gs and integration of 
the differential (|39|) one gets: 

r v 2 N- 

W^gaspfey*- (41) 

Here v = v(j) = 7 - 7^ and v in = v(j [n ). And p(v Pj ) = pj , j = 1, 2. 

The solution of (|4T|) is given in terms of the elliptic <r and ( functions: 



, / n 1 - V 7 ^ |«7) tA . 

^ (7) o arctan / J ( 43 ) 

^ 4\Er. 7=777. /l _ 1, 2 



can be easily integrated and the solution for < and D > is given by 

where is defined in ([36]) . 

4.5 Solution of the t— equation 

Like for the coordinates tp and i/j the Eq. ([9]) can be simplified by using ((5]) and by performing the substitution 
u = ^ (4y — y) (for a = 6 only the ix-part survives) 

^ = ^=fe^o + E^) ■ (^) 

where the partial fractions decomposition procedure was applied. Here P2 and p\ are defined as before in the <p and 
ip equations. T^, i = 0, 1, 2 are constants which arise from the partial fractions decomposition. These depend on the 
parameters of the metric and the test particle. 
The integration of ([44]) yields: 

t-t in = t u , (45) 

where t\ n is the initial value of t. 



Substitution of y = p(u) with p'(v) = ^/4p 3 (v) — g2p(v) — #3 in (|44|) yields: 



7} 



0) " PK>J 



. (46) 



Here w = ^(7) = 7 - 7^ and v in = v(y m ), and p{v Pj ) = pj , j = 1, 2. 

The solution of (jl"6")) is given in terms of the elliptic <r and £ functions: 

,„ = £ («„,„) - «,,„ + To( „ - „,„) + ± _i_ U p>)( „ _ + h ^Sj - ,» , (47) 



5 The orbits 



In this section we show examples of the types of orbits discussed in the previous sections. The orbits are plotted in 
the cartesian coordinates (X, Y, Z, W) defined in [31] as 

X = yj u + a 2 sin d cos <p , Y = \j u + a 2 sin d sin , (48) 
Z = \/u + a 2 cos $ cos -0 , W = y/u~+~a? cos $ sin ^ , (49) 

where r G (0, oo) , $ G (0, ^) , G (0, 27r) , ^ G (0, 27r). Since we cannot visualize a 4D picture we present 3D orbits 
for e.g. the coordinates (X, Y, Z) in the Figj6j Every 3D orbit is bounded by two cones with the opening angles 
defined by the roots of the polynomial ([14|) . The spheres in the plots identify the outer and inner horizon. 

Because of the bad choice of the coordinates the geodesies diverge at the horizons. This is associated with the 
directional change of the orbits in Fig. E ^a)|(e)|(T)| and |(h)| It can be already seen in the differential equations ([7j) 
and (|HJ) because of the quadratic polynomial A which defines the horizons. This behaviour is known from the 4 
dimensional Kerr(- de Sitter) space-time [9j[in|[23]. In the next section [6] we proof that the frame-dragging effect is 
non zero in the 5-dimensional Myers-Perry space-time. This can be also seen in the MBO in the Figj ^a)| and |(h)| 
BO in the plots |(c)| and |(d)| as well as in the escape orbit |(g)| Here the orbit in the plot |(c)| has the energy value 
close to the energy of an unstable circular orbit in the potential well of the potential |4(f)| 

Another interesting feature of the motion around rotating black holes can be at best seen in the pictures |6(e 



and |(h) | (and also in the plots |6(a)] and [6(f)] but feebly seen). When the angular momentum (or both of them) of 
a test particle is negative in contrast to the positive angular momentum of the black hole, a test particle will be 
dragged in the direction of the black hole's rotation: as soon as it approaches the ergosphere it starts to corotate 
with the black hole. In the pictures |6(e) and |(h)| this behaviour correponds to the kinks shortly before the horizon 



(i.e. at the ergosphere which is not plotted here). That does not happen if one of the angular momenta has the 
same sign as the rotation parameter of the black hole (see e.g. picture |6(f)| ) . 

In the Figdwe show a bound and an escape orbit for massless particles illustrating similar properties as described 
above for the massive test particles. 



6 The observables 

For observations one can define observables which do not depend on the chosen coordinate frame. These are e.g. 
the perihelion shift for bound orbits, the light deflection for escape orbits, the deflection angle for fly by orbits, or 
the Lense-Thirring effect which is present in axially symmetric space-times. We follow along the lines of [32 , 33] to 
calculate the observables. 

Consider the perihelion shift for a bound orbit BO or many world bound orbits MBO. 

The it-motion is periodic with period 



du f e2 dy 



n VU J ei JP^y) ' 

2-7T 

where e\ and e2 are the zeros of P^{y) related to ii m in and u max . The corresponding orbital frequency is — . 

2tt 

The polar frequency is given by — , where the polar period of the ^-motion is defined as 



dti n p max df 

"* = 2 / 7^ = " 2 



for the angle tf min < i? < tf max where tf min ,^max e (0,tt). 

Next we calculate the secular accumulation rate of the time t: 



where t u is defined in (|47|) and j ei corresponds to the root e^, i = 1,2 
The secular accumulation rates of the angles cp and ip are given by: 



(50) 



fMnax fit 

I ~w (51) 



7e 2 
7ex 



(53) 



(a) MBO. r = l,a = 0.2, K = 3,E = 
1.14756215, <£ = 0.2, ^ = 0.2. 



(b) EO. r = l,a = 0.2, K = 3, £ = 
1.14756215, $ = 0.2, ^ = 0.2. 




(c) BO. r = l,a = 0.4, AT = 10, £ = (d) BO. r = l,a = 0.4, K = 10, E = (e) TEO. r = 1, a = 0.4, K = 10, E = 

8.16516573, 3> = -4, ^ = -0.5. 8.126604544, $ = -4, ^ = -0.5. 8.126604544, $ = -4, ^ = -0.5. 




(f) TEO. r = l,a = 0.2, X = (g) EO. r = I, a = 0.4, X = 5, £ = (h) MBO. r = l,a = 0.4, if = 5, £ = 
3, £7 = 1.18247681,$ = 1.625363954, 3> = -1, ^ = -1. 1.625363954, $ = -1, ^ = -1. 

-0.3, ^ = 0.9. 



Figure 6: Orbits for 6 = 1 corresponding to the possible types specified in the Tabled MBO, BO, EO and TEO. In the picture [(e)] we omit the 
upper cone for better view. 



and 



7e 2 
7ei 



(54) 



(a) BO. r = l,a = 0.4, K = 10, E = (b) EO. r = l,a = 0.4, AT = 5, £ = (c) EO. r = l,a = 0.4, AT = 1, £ = 

8.07660455, $ = -4, ^ = -0.5. 1.521633952, 3> = -1, ^ = -1. 0.52, $ = 0.01, ^ = -0.01. 



Figure 7: Orbits for 5 = co rresponding to the possible types specified in the potential plots [5] E.g., the depicted bound orbit in |(a)| relates to 
the potential |5(e)] Spheres have the radia of the horizons. 



The orbital frequences £l u , Q$ and f^, are then given by: 

o 2lTl o 27Tl o y * o Y * m\ 

UJ U I Co># I 1 1 

The perihelion shift and the Lense-Thirring effects are defined as differences between these orbital frequences 

^ _ Q _ Q -l( Y - —\ 

^Perihel ~ J V Si « ~ p 1 ^ ^ I 



1 2tt 



^Perihel ~ ~ r ^ 



r ^ w„ 



(56) 



^Lense-Thirring ~ "V # T \ 



^Lense-Thirring — ^ — p ( ^ 



1 (, r 2tt\ 



- 



It seems to be not possible to measure a perihelion shift for a particle moving on an orbit located inside the 
inner horizon or for a many world bound orbit where only one part of the orbit can be seen and other parts are 
believed to happen in other universes. A deflection of light or massive particles by the gravitational field of the 
black hole seems to be more likely observable. This can be calculated in an analogous way (here in the ^-motion 
the upper limit is infinity). 

7 Conclusions and Outlook 

In this paper we presented the analytical solution of the equations of motion in the 5D Myers-Perry space-time with 
equal rotation parameters. We integrated the encountered differentials of the first and third kind in terms of the 
Weierstrass's elliptic, zeta- and sigma-functions. We studied the properties of test particle motion defined by the d- 
and u- polynomials and completely characterized the possible types of motion. As examples we plotted the orbits 
of massive test particles for different values of the characterizing parameters. 

The next step is to investigate the motion around a Myers-Perry black hole with non-equal rotation parame- 
ters [8]. It will elucidate the influence of the rotation of the black object on the motion of test particles further. 
It would be also interesting to go to higher dimensions, to add NUT parameters and to include a cosmological 
constant. This will increase the order of the polynomials and the solution will be expressed in terms of hyperelliptic 
functions. 
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